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Description 

BACKGROUND OF THE INVENTION 

(1) Field of the Invention 

[0001] The present invention relates to a sound envi- 
ronment simulator which enables one to experience a 
computer simulation of sound performance when mod- 
eling acoustic affects in a listening room, a concert hall 
or the like, and to a method of analyzing spatial acoustic 
characteristics of a sound space. 

(2) Description of the Related Art 

[0002] Rapid advancement in the field of audio-visual 
technology makes it possible to supply high-quality mu- 
sic and images to individual homes. In particular, what 
has been attracting people's interest is a technique, a 
so-called "home theater", by which one can enjoy a 
movie or music at home as if he were in a movie theater 
or a concert hall. To enjoy the home theater, it is desir- 
able to have a room where acoustic affects and sound 
insulation effects are well taken into account. However, 
in most of the cases, people enjoy movies and music in 
their living rooms where a variety of furniture and house- 
hold appliances are installed. Thus, nowadays, not only 
the appearance and the feeling of being comfortable, 
but also the acoustic affects and insulation effects are 
important elements when designing for living rooms. 
[0003] Conventionally, data related to acoustic char- 
acteristics are evaluated for modeling acoustic affects: 
they are computed by either experimental or analytical 
method and then fed back to design specifications. 
[0004] In the experimental method, a model room is 
constructed, and an impulse response is measured di- 
rectly to compute the acoustic characteristics data of the 
model room. However, the costs both for building mate- 
rials and labor do not allow its application to individual 
house construction. 

[0005] Whereas in the analytical method, the acoustic, 
characteristics data are computed by a computer as it 
simulates a three dimensional sound field. This method 
has been widely used in recent years. Because the re- 
sulting data reliably correspond to various input condi- 
tions, become a reliable source for further analysis, and 
facilitate further processing and feedback to the design 
specifications. 

[0006] Data necessary for simulating the three dimen- 
sional sound field are computed either by a numerical 
analysis or a geometric analysis based on theory of 
acoustics: sound waves are taken into account in the 
former, whereas they are not in the latter. 
[0007] A finite element method(FEM) and a boundary 
element method(BEM) are typically used in the numer- 
ical analysis. In both the methods, the computer com- 
putes acoustic intrinsic modes or pressure distribution 
of a sound by solving a dominance equation of a sta- 



tionary sound field, or namely Helmholtz's equation. 
Taking the sound waves into account makes it possible 
to analyze diffraction and interference of the sound; 
however, it does not allow the computation of a non-sta- 

5 tionary value such as an echotime pattern. In particular, 
to compute the echotime pattern of a high frequency 
sound, a myriad number of lattices are involved, so that 
not only it takes relatively long, but also the amount 
thereof exceeds a memory capacity. For this reason, the 

io numerical analysis has not been applied to practical ap- 
plications. 

[0008] On the other hand, in the geometric analysis 
based on the theory of acoustics, methods using sound 
rays or virtual images are typically used on the premise 

is that the sound rays are reflected geometrically from 
walls. These methods are practical; for the memory ca- 
pacity does not limit the amount of the computation. 
However, this analysis does not guarantee accuracy in 
ail cases. For example, the resulting data may not be 

20 accurate when simulating the sound field with a low fre- 
quency sound in a room where a wave length is longer 
than its size. Because given these circumstances, the 
sound rays fail to reflect in accordance with the premise. 
[0009] In addition, a reverberation time, which is a crit- 

25 jcal element to compute the acoustic characteristics da- 
ta, is approximated by Sabine's or Eyring's equation to 
minimize the amount of computation for the sound ray 
tracing. Such an approximation causes an computation 
error when this analysis is applied to a space where a 

30 variety of household appliances are to be installed or 
whose shape is complex. 

[0010] Another type of sound environment simulator 
has been developed in recent years. This type of simu- 
lator simulates a sound field by reproducing a sound in 

35 accordance with the acoustic characteristics data, so 
that one can hear the reproduced sound. In a broad 
sense, the sound environment simulator includes a 
model room where one can listen to a sound reproduced 
over speakers installed therein. However, this type of 

40 simulator is not applicable for individual house construc- 
tion. 

[0011] Accordingly, a sound environment simulator, 
disclosed in U.S. patent application No. 4,731,848, was 
proposed. This simulator is designed to reproduce a 

45 sound over a headphone with reverberation effects. 
More precisely, an input audio signal is branched and 
the branched signal is inputted into a delay circuit. Then, 
the delay circuit delays the output of the branched signal 
for latency comparable to the time difference of the re- 

50 fleeted sound to produce a reverberant signal. Subse- 
quently, the original input audio signal and reverberant 
signal are composited to generate an output signal 
which is transduced into a sound. However, the sound 
reproduced in this way may give an impression to a lis- 

55 tener that it does not sound realistic; for the sounds of 
higher order reflections (i.e. the sounds reflected for a 
number of times) are discarded. 
[0012] Moreover, this type of simulator does not re- 
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produce the sound immediately. This is because index 
data such as locations of a virtual sound source and the 
listener are inputted manually each time they move. Fur- 
thermore, the shape of a space selected as a model en- 
vironment is limited. 

[0013] US-A-4 928 250 discloses a system for deriv- 
ing radiation images of an environment with an algo- 
rithm using form factors for setting analysis conditions. 
To find the form factors, a hemi-cube is constructed 
around the surface with grid cells defined for all faces 
on the hemi-cube. The identity of the polygons in the 
environment are stored to speed up the calculations. 
While US-A-4 928 250 mentions the possibility of using 
its invention for sound, US-A-4 928 250 discloses only 
radiation images of light in details and does not disclose 
how to apply it to sound. Since light radiation does not 
involve noticeable time delay, US-A-4 928 250 does not 
disclose or suggest use of an echotime pattern which 
represents time delays of radiation reverberation. 
[0014] JP-A-11 041 321 discloses an audible type 
sound field simulator which uses a signal into a rever- 
beration adding device 14 branched from the sound 
source music 9 to produce reverberation sound 1 5. The 
reverberation sound is added to music data 11, 12 to 
obtain music data 16, 17. This device uses an echo time 
pattern obtained by acoustic wave simulation by a virtual 
image method, and does not use a geometric method 
based on the theory of acoustics. 
[0015] JP-A-30 038 695 discloses an audible in-room 
sound field simulator comprising an initial reflected 
sound generation part, an initial scattered sound gener- 
ation part and a reverberation generation part. The in itial 
reflect sound I, initial scattered sound S, reverberations 
R are generated separately by the respective parts of 
the simulator. This simulator does not use an echo time 
pattern to compute reverberation sound or initial scat- 
tered sound. 

SUMMARY OF THE INVENTION 

[0016] Accordingly, the present invention has an ob- 
ject to provide a sound environment simulator which en- 
ables a listener to experience a simulation of sound en- 
vironment of a space selected as a model environment 
without constructing the model space. 
[0017] According to an aspect of the present inven- 
tion, there is provided a sound environment simulator 
which comprises echotime pattern computing means for 
computing an echotime pattern of a space; impulse re- 
sponse computing means for computing an impulse re- 
sponse using said echotime pattern; reproduced sound 
generating means for generating reproduced sound da- 
ta using said impulse response; and output means for 
outputting a sound by transducing said reproduced 
sound data into the sound. The present simulator is 
characterised in that said echotime pattern computing 
means is adapted to compute said echotime pattern by 
a geometric method based on the theory of acoustics 



with input data as to conditions of said space; said im- 
pulse response computing means is adapted to com- 
pute a response for each frequency band from said 
echotime pattern and to compute a composite response 

5 to generate said impulse response; and said repro- 
duced sound generating means is adapted to generate 
the reproduced sound data by convolution of said im- 
pulse response on non-reverberate sound source data. 
[0018] The simulator may further comprise analysis 

10 data setting means for setting analysis data by receiving 
analysis condition data including a shape and physical 
data of walls of said space, dividing said walls into a set 
of certain sections, and setting a virtual sound source 
and a virtual sound receiving point at respective arbi- 

1$ trary positions within said space. Said echotime pattern 
means may be adapted to receive said setting analysis 
data and include time series data computing means for 
computing, based cn said setting analysis data, time se- 
ries data related to arrival volume of a sound emanated 

20 from said virtual sound source to said virtual sound re- 
ceiving point; and data transducing means for transduc- 
ing said time series data related to the arrival volume 
into said echotime pattern by interpolating certain data. 
[001 9] The simulator may further comprise dry source 

25 input means for inputting said non-reverberate sound 
source data; and head transfer function storing means 
for storing data related to a head transfer function; 
wherein said reproduced sound generating means is 
adapted to generate the reproduced sound data by con- 

30 voluting said impulse response received from said im- 
pulse response computing means and said head trans- 
fer function received from said head transfer function 
storing means on said non-reverberate sound scurce 
data received from said dry source input means. 

35 [0020] The simulator may further comprise location 
data detecting means for detecting data related to at 
least one of a listener's head location and a listener's 
head direction; wherein said reproduced sound gener- 
ating means is adapted to extract said impulse response 

40 and said head transfer function corresponding to said 
detected location data to convolute them on said non- 
reverberate sound source data. 
[0021] The location data detecting means may be 
adapted to detect time varying data. The output means 

45 may comprise a headphone. 

[0022] The simulator may also comprise head trans- 
fer function storing means for storing data related to a 
head transfer function; and composite response gener- 
ating means for computing a composite response for 

so each direction of a listener's head by convolution of said 
impulse response and said head transfer function, 
wherein said reproduced sound generating means is 
adapted to generate reproduced sound data by convo- 
lution of said composite response on the non-reverber- 

S5 ate sound source data. 

[0023] The time series data computing means may be 
adapted to compute form factors to compute the arrival 
volume of the sound to the sound receiving point by us- 
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ing said form factors and a sound absorption coefficient 
of each section included in said physical data of the 
walls. 

[0024] According to another aspect of the present in- 
vention, there is provided a methcd of analysing a sound s 
space, comprising the steps of computing an echotime 
pattern of a space; computing an impulse response us- 
ing said echotime pattern; generating reproduced sound 
data related to a reproduced sound based on said im- 
pulse response; and outputting a sound by transducing 
said reproduced sound data in the sound. The method 
is characterised by the step of detecting data related to 
a listener's location; wherein the step of computing an 
echotime pattern computes said echotime pattern for a 
virtual sound source placed arbitrarily in the space by a 
geometric method based on the theory of acoustics: the 
step of computing an impulse response computes said 
impulse response by convoluting said echotime pattern 
for each frequency band by a bandpass filter to compute 
a response for each frequency band, and by computing 
a composite response to compute said impulse re- 
sponse; and the step of generating reproduced sound 
data generates said reproduced sound data by receiv- 
ing a dry source of the sound source, and by convoluting 
said impulse response and a head transfer function cor- 
responding to the detected listener's location on said dry 
source of the sound source. 

[0025] In the step of detecting data related to the lis- 
tener's location, data related to dynamic changes of the 
listener's position and the direction of his head may be 
detected as the location data. 

[0026] In the step of detecting data related to the lis- 
tener's location, data related to directional changes of 
the listener's head may be detected as the location data. 
[0027] The method may further comprise the step of 
dividing walls of the space into a set of sections; setting 
analysis conditions including the shape of the space, 
and physical data of the walls, the set of sections, a vir- 
tual sound source and a virtual sound receiving point; 
wherein the step of calculating said echotime pattern 
may include the steps of computing time series data re- 
lated to arrival volume of a sound emanated from said 
sound source to said sound receiving point for each in- 
cident direction and frequency band; and computing 
said echotime pattern by interpolating certain data to 
said time series data. 

[0028] The step of computing time series data may 
include the steps of computing data related to the arrival 
volume of the sound emanated directly from said sound 
source to said sound receiving point; computing data re- 
lated to the arrival volume of the sound arriving at said 
sound receiving point by reflecting for up to a predeter- 
mined number of times by using a sound absorption co- 
efficient and form factors of each section included in said 
physical data of the walls; computing data related to the 
arrival volume of the sound arriving at said sound re- 
ceiving point by reflecting for more than the predeter- 
mined number of times by using the sound absorption 



coefficient and the form factors of each section included 
in said physical data of the walls; and computing time 
series data by dividing said arrival volume data into a 
set of groups by a certain time interval. 
[0029] In the step of computing data related to the ar- 
rival volume of the sound arriving at said sound receiv- 
ing point by reflecting for more than a predetermined 
number of times, a new virtual sound source may be set 
by integrating said arrival volume data of each group, 
and the arrival volume of sounds emanated from said 
new virtual sound source to any of the sections may be 
repeatedly computed. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0030] These and other objects, advantages and fea- 
tures of the invention will become apparent from the fol- 
lowing description thereof taken in conjugation with the 
accompanying drawings which illustrate specific em- 
bodiments of the invention. In the drawings: 

Fig. 1 is a block diagram of the sound environment 
simulator of the present invention; 
Fig. 2 is a block diagram of the sound environment 
simulator in accordance with the first embodiment; 
Fig. 3 is a flowchart detailing the operation of the 
sound environment analyzing unit 1 ; 
Fig. 4 is a flowchart detailing the sub-routine of the 
step for computing the echotime pattern in Fig. 3; 
Fig. 5 is an illustration depicting the model space 
used by the echotime pattern computing unit 15; 
Fig. 6 is another illustration depicting the model 
space used by the echotime pattern computing unit 
15; 

Fig. 7 (a) is a graph of the time series data of the 
actual arrival volume at one of the divided section; 
Fig. 7(b) is another graph of the time series data of 
the arrival volume divided into a set of groups at one 
of the divided section; 

Fig. 8(a) is a graph of the time series data of the 
actual arrival volume at the sound receiving point; 
Fig. 8(b) is a graph of the time series data of the 
arrival volume divided into a set of groups at the 
sound receiving point; 

Fig. 8(c) is a graph of the time series data of the 
arrival volume interpolated with a certain data 
(echotime patter) at the sound receiving point; 
Fig. 9 is a block diagram detailing the construction 
of the sound field reproducing unit 2 and output unit 
3; 

Fig. 10 is a block diagram detailing the construction 

of the impulse response computing unit 21; 

Fig. 11 is a flowchart detailing the operation of the 

impulse response computing unit 21; 

Fig. 12 is a flowchart detailing the operation of the 

sound field reproducing unit 2 in Fig. 9: 

Fig. 1 3 is a block diagram of the sound environment 

simulator in accordance with the second embodi- 



es 



20 



25 



30 



35 



40 



45 



SO 



4 



7 



EP 0 593 228 B1 



8 



ment; and 

Fig. 1 4 is a block diagram of the sound environment 
simulator in accordance with the third embodiment. 

DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 

[First Embodiment] 

[0031] A sound environment simulator in accordance 
with the first embodiment reproduces a realistic, stere- 
ophonic sound by computing data related to acoustic 
characteristics in a virtual space created as a model en- 
vironment, so that a listener can experience a simulation 
of a sound field. 

[0032] The sound environment simulator, as is shown 
in Fig. 1 , comprises a sound environment analyzing unit 

I, a sound field reproducing unit 2, and an output unit 
3. The sound environment analyzing unit 1 analyzes the 
acoustic characteristics in the virtual space by means of 
computer simulation to compute a echotime pattern. 
The sound field reproducing unit 2 computes data relat- 
ed to an impulse response from the echotime pattern to 
transduce non-reverberant portion of a desired sound 
into a consecutive sound signal on real time. The output 
unit 3 includes a headphone or the VR equipment and 
transduces the consecutive sound signal into a sound. 
[0033] The above-described sound environment sim- 
ulator is designed to simulate changes in the sound oc- 
curring as the listener stands up, sits down, walks 
around in a room, or moves his head. More precisely, 
the sound environment analyzing unit 1 computes the 
echotime patterns for each location, incident direction, 
and frequency of a sound source. The sound field re- 
producing unit 2 computes the impulse response from 
the echotime patterns, and adds further processing to 
produce a stereophonic sound signal, which is trans- 
duced into a sound by the output unit 3. 

[0034] Fig. 2 provides a detailed construction of the 
sound environment analyzing unit 1. it comprises a 
room data input unit 11, a sound source data input unit 
1 2, a sound receiving point data input unit 1 3, a surface 
dividing unit 14, an echotime pattern computing unit 15 
including a time series data computing unit 1 6 and a da- 
ta transducing unit 17, and a memory 18. 
[0035] The sound environment analyzing unit 1 com- 
putes the echotime patterns for each location in the vir- 
tual space, each incident direction of the sounds inci- 
dent to the sound receiving point, and each sound 
source including a plurality of predetermined frequency 
bands. 

[0036] The flowchart in Fig. 3 details the above-de- 
scribed operation. 

[0037] In Step 100, data related to the virtual space 
such as the shape, size and a sound absorption coeffi- 
cient of walls are inputted into the room data input unit 

II. Subsequently, data related to the location of the 
sound receiving point(location of the listener's head), 



and data related to incident direction of the sound to the 
sound receiving point are inputted at the sound receiv- 
ing point data input unit 13 in Steps 110 and 120, re- 
spectively. Further, data related to the sound source 

5 such as the frequency and volume thereof are inputted 
at the sound source data input unit 12 in Step 130. 
[0038] Then, in Step 140, the surface dividing unit 14 
divides solid surfaces, or namely the walls, of the virtual 
space into a set of sections, and numbers them serially. 

10 Each section is a rectangle whose side is shorter than 
the wave length of the sounds emanated from the sound 
source. It also sets the sound absorption coefficient for 
each section from the sound absorption coefficient in- 
putted in Step 100, and the sound receiving point. The 

1$ sound receiving point takes a form of, for example, a 
dodecahedron to approximate the listener's head. 
[0039] An example of the virtual space created in this 
way is illustrated as a model space 10 in Fig. 5. In the 
drawing, one wall is divided into approximately 400 sec- 

20 tions, and the sound source is set at 500 Hz, whereby 
the short side of each section is set to 45cm. The sound 
absorption coefficient a is added to each section; the 
sound absorption coefficient ai is added to the section 
i, and aj to the section j. The frequency of the sound 

25 source is set per 1/8 octave band; forthe sound absorp- 
tion coefficient data of the walls are practically computed 
per octave band. 

[0040] Once the model space 1 0 is created, the echo- 
time pattern calculator 15 computes the echotime pat- 
30 tern in Step 200, which is detailed by the flowchart in 
Fig. 4. 

[0041] In Step 210, the time series data computing 
unit 16 computes form factors. These form factors are 
calculated by a radiant heat ray tracing technique, which 

35 is often used in analyzing radiation in the thermal envi- 
ronment or lighting. Japanese laid-open patent applica- 
tion No. H5-52383, Hattori et al teaches the calculation 
method in detail. In this step, the form factors are calcu- 
lated according to this calculation method with respect 

40 to the relations between: the sound source and each 
section; one section and another section; and each sec- 
tion and the sound receiving point, which are respec- 
tively referred to as Fsc, Fee and Fcr. 
[0042] The calculation of reflected sound with the 

45 form factors further adds reflection components having 
directivity and those of mirror reflection besides perfect 
reflection components to form factors as is with the light- 
ing. Thus, more realistic reflections can be simulated by 
using this calculation. Also, this calculation is virtually a 

so must when only an initial reflection by the mirror reflec- 
tion at a relatively large divided section of a huge space 
such as a concert hall is roughly estimated, or when the 
reflection has the directivity due to the inclined walls. 
[0043] More precisely, each form factor is calculated 

55 in the following way: 

1) Form Factor Fsc (Sound Source -> Section) 
Suppose an arbitrary positioned virtual sound 
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source S evenly emanates a plurality of sound rays 
(e.g. 100,000 rays) in all directions. Then, specific 
sections at which each sound ray arrives are detect- 
ed. For example, with the model space 10 in Fig. 5, 
a sound ray SRsk arrives at a section k, and sound 
rays SRsi 1 and SRsi 2 arrive at a section i. Accord- 
ingly, the time series data computing unit 16 com- 
putes the arrival volume P1 at each section to com- 
pute P1 : Ps, PVs ratio to the volume of the sound 
source Ps, which is referred to as the form factor 
Fsc. In this way, the form factors Fsc's (or all the 
sections are calculated. 

2) Form Factor Fee (Section Section) 

For factor Fee is calculated in the same manner. 
In this case, the arrival volume to one of the sections 
is considered to be a new sound source. 

3) Form Factor Fcr (Section -> Sound Receiving 
Point) 

Form factor Fcr is calculated in the same man- 
ner. 

[0044] Once these form factors are calculated, the 
time series data computing unit 1 6 computes time series 
data related to the arrival volume. Note that the sound 
receiving point R receives two types of sounds : one is 
a direct sound from the sound source, and the other is 
a reflected sound from one of the sections. Thus, the 
time series data computing unit 16 computes the time 
series data of the direct sound first, thence, those of the 
reflected sound. For example, with the model space 10, 
the sound emanated along with the sound ray SRsr is 
the direct sound, and the one with the sound ray SRsk 
is the indirect sound. 

[0045] Back to the flowchart, the time series data 
computing unit 16 computes the arrival volume of the 
direct sound in Step 220. The arrival volume denoted as 
Pr is computed with the distance from the sound source 
to the sound receiving point R according to the inverse 
square law, while at the same time, the arrival time de- 
noted as TR is computed by dividing the distance by the 
speed of sound. 

[0046] Next in the Step 230, the time series data com- 
puting unit 1 6 computes the arrival volume of the reflect- 
ed sound. The amount of this computation increases a 
geometric series-fold, because the number of reflected 
sound rays are multiplied by the number of the sections 
each time they reflect. Thus, a memory must have quite 
a large capacity to store these myriad computation re- 
sults. To avoid this problem, the volume of the reflected 
sound is computed in two ways depending on the 
number of reflections, or largeness of the orders. In this 
embodiment, the reflections up to the third order are re- 
ferred to as lower order reflections, and the reflections 
beyond the third order reflections as higher order reflec- 
tions. 

[0047] To begin with, a method for computing the ar- 
rival volume of the sound of the lower order refections 
to the sound receiving point R will be explained with re- 



ferring to Fig. 5. 

(The First Order Reflections) 

s [0048] Let the sound rays SRsk and SRkr be exam- 
ples. The sound along with the sound rays are emanat- 
ed from the sound source S and incident to the sound 
receiving point R by being reflected once at the section 
k. The volume thereof denoted as Pr 2 is computed with 

io the form factors Fsc-, and F^r, and absorption coeffi- 
cient cxk according to : 



25 



35 



45 



Pr 2 = Fsc 1 x Fc 1 r x (1 - cck) x Ps 



15 



[0049] Subsequently, like with the direct sound, the 
distances between the sound source S and the section 
k, the section k to the sound receiving point R are added 
and divided by the speed of sound to compute the arrival 
20 time T 2 . 

[0050] The above calculation is applied to all of the 
sounds of the first order reflections. 



(The Second Order Reflections) 



[0051] Let the sound rays SRsi 2 , SRij and SRjr be ex- 
amples. The sound along with these sound rays are re- 
flected at the sections i and j before it arrives at the 
sound receiving point R. The arrival volume denoted as 
30 Pr 2 is computed with the form factors Fsc, Fccg and Fc 2 r, 
and absorption coefficients ai and ak according to : 



Pr 2 = Fsc x Fcc 2 x Fc 2 r x (1 - ai) x (1 - ak) x Ps 



[0052] Also, the distances between the sound source 
S and the section i, the sections i and the section j, and 
the section j to the sound receiving point R are added 
and divided by the speed of sound to compute the arrival 
40 time T 3 . 

[0053] The above calculation is applied to all of the 
sounds of the second order reflections. 



(The Third Order Reflections) 

[0054] The volume of the sound of all of the third order 
reflections is computed in the same manner described 
as above. 



so (The Higher Order Reflections) 

[0055] Next, a method for computing the arrival vol- 
ume of the sound of the higher order refections to the 
sound receiving point R will be explained with referring 
55 to Fig. 6. 

[0056] Let the section n be an example, then the time 
series data of the third order reflections, SRjn and SRmn 
as is shown in Fig. 7(a), have been computed at the end 
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of the computation for the third order reflections. 
[0057] These time series data are divided into a set 
of groups by a certain time interval At, for example, 0.1 
Ms. The volume of each group is added to find a sum of 
the time series data within each group, and the resulting 
data, referred to as the time series data of the sum, is 
shown in Fig. 7(b). Then, the section having the time 
series data of the sum is positioned as a new virtual 
sound source S'. In practice, at the new sound source 
S', N sound rays arriving at the concerned section during 
At are composited into a single sound ray, so that it will 
have a volume equal to the sum of the volume of the N 
sound rays'. 

[0058] Back to the flowchart again, in Step 250, as is 
with the sound of the lower order reflections, the arrival 
volume along with the sound rays to either the sound 
receiving point R or another section are computed. This 
computation is repeated until the arrival volume of all 
the sections becomes 1/10000 of the entire emanated 
volume Ps or less. 

[0059] During the above-described computation, the 
sound receiving point R receives a large amount of the 
direct and indirect sounds as is shown in the graph in 
Fig. 8(a). The arrival sound rays are divided into a set 
of groups as is with the higher order reflections, for ex- 
ample by a certain time interval, then the number thereof 
is counted, and the arrival volume of each sound ray are 
summed up if they are plural. 

[0060] When the arrival volume of all the sections be- 
comes lower than the predetermined level in Step 250, 
the further computation is not proceeded. Instead, such 
a petty arrival amount are summed up in each section, 
and these sums are further added to the time series data 
in Step 260. 

[0061] As a result, the time series data of the arrival 
volume at the sound receiving point R, as is shown in 
Fig. 8(b), is computed. 

[0062] Accordingly, the data transducing unit 1 7 com- 
putes a mean arrival volume in Step 270 to interpolates 
data to generate the echotime pattern from the time se- 
ries data of the sum in the following step 280. 
[0063] More precisely, the data transducing unit 17 
computes a mean arrival volume data as is shown in 
Fig. 8(c) from the time series data of the sum. To com- 
pute the mean arrival volume data, the sum of the arrival 
volumes at the receiving point R within a certain group 
is divided by the number of the arrival sound rays. As- 
suming that a composite sound ray SRr consisting of N 
(=3) sound rays arrives to the sound receiving point R 
with a arrival volume Pm during the third group(0.2-0.3 
Ms), then the mean volume data SRa is computed ac- 
cording to : 

SRa = Pm/N 

[0064] Subsequently, to interpolate data to the mean 
volume data SRa thus calculated, the data transducing 



unit 17 applies the mean volume data of each group to 
the intermediate point of the corresponding group. 
Then, as is shown in Fig. 8(c), the tops of each bar in 
the graph representing the mean volume data of the 
s groups are finked, so that an approximate line will be 
formed. By using this approximate line, the data con- 
verting unit 17 interpolates as many new lines as the 
actual arrival N sound rays in the corresponding group. 
As a result of this data interpolation, the echotime pat- 
10 tern of an approximate line as is shown in Fig. 8(c) is 
formed. This interpolation is not limited to the approxi- 
mation by a line; a curve may be used instead of a line, 
or an addition of a weight corresponding to the arrival 
time to each sound ray also will do. 

75 [0065] The optimal effects are obtained when At is set 
to a corresponding value to a sampling frequency of the 
audio equipments such as compact disc player; for ex- 
ample, the optimal At for a sampling frequency of 44.1 
KHz is 22.7 Ms. 

20 [0066] The echotime pattern thus computed is stored 
into the memory 18, and read out for further processing 
at the sound field reproducing unit 2. 
[0067] Fig. 9 provides a detailed construction of these 
two components. The sound field reproducing unit 2 

25 comprises an impulse response computing unit 21 for 
reading out the echotime pattern from the memory 18 
to compute an impulse response, a location data input 
unit 22 for reading out the data of the location of the 
listener within the model space; the location data are 

30 detected by a location data detecting unit 31 such as the 
VR equipment in the output unit 3, a head transfer func- 
tion storing unit 23 for storing a predetermined head 
transfer function, a dry source input unit 24 for inputting 
data related to a dry source of a sound(the sound re- 

35 corded in an anechoic room) to be reproduced, a repro- 
duced sound generating unit 25 for generating a con- 
secutive sound signal by convoluting data related to the 
impulse response and head transfer function on the data 
related to the dry source with a plurality of digital signal 

40 processors(DSPs). 

[0068] As is shown in Fig. 10, the impulse response 
computing unit 21 includes a echotime pattern input unit 
41 for inputting data related to the echotime pattern for 
each location in the model space, each incident direc- 

45 tion, and each frequency band, a convolution computing 
unit 42 for convoluting the echotime pattern of each fre- 
quency band by a bandpass filter, an impulse response 
computing unit 43 for computing the impulse response 
by compositing responses from all the frequency bands, 

50 an impulse response output unit 44 for outputting the 
impulse response for each location and incident direc- 
tion, and a memory 45 for storing data related to the im- 
pulse response. 

[0069] The output unit 3 includes a sound output in- 
55 strument 32 such as a headphone, over which a sound 
based on the consecutive signal from the reproduced 
sound generating unit 25 is reproduced, so that the lis- 
tener can hear it. 
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[0070] In practice, a RB2(Reality Built for 2, a trade 
name) system or an Eyephone system of VPL Research 
Inc, and POLHEMUS are used for the location data de- 
tecting unit 31. In the Eyephone system, a sensor placed 
on top of the listener's head detects his location and his 
head's direction to output the related data to the sound 
field reproducing unit 2 on real time. 
[0071] The above described sound field reproducing 
unit 2 and output unit 3 operate in the following way. 
[0072] The impulse response computing unit 21 com- 
putes the impulse response for each location and inci- 
dent direction (listener's head direction) from the echo- 
time pattern. 

[0073] The flowchart in Fig. 11 details how the impulse 
response is generated from the echotime pattern. 
[0074] in Step 300, the impulse response computing 
unit 21 reads out the echotime pattern sequentially from 
the memory 1 8, and in Step 320, it computes a response 
for each frequency band by convoluting each frequency 
by a bandpass filter. Consequently, the impulse data re- 
lated to all the locations in the model space and incident 
direction are computed and stored in the memory 45. 
[0075] Once the impulse data are computed, the 
sound field reproducing unit 2 shifts to a sound repro- 
ducing process, which is detailed by the flowchart in Fig. 
12. 

[0076] In Step 400, the sound field reproducing unit 2 
reads out the location data from the location data input 
unit 22 on real time. Then, the impulse response input 
unit 41 reads out the impulse response data corre- 
sponding to the location data from the memory 45 in 
Step 41 0. Further, it reads out the head transfer function 
data from the head transfer function storing unit 43 cor- 
responding to the location data as to the listener's head 
direction in Step 420. As previously mentioned, the head 
transfer function is a parameter that represents the 
change in an incident sound in accordance with the 
shape of the listener's head and a location of his ears, 
and stored in the head transfer function storing unit 43 
in advance. 

[0077] Subsequently, in Step 430, the reproduced 
sound generating unit 25 receives the dry source of the 
sound from the dry source input unit 24, and convolutes 
the impulse response data on head transfer function da- 
ta with the DSPs to generate a consecutive sound sig- 
nal, which is outputted to the sound output unit 52 on 
real time and transduced into a sound in Step 440. 
[0078] As a result, the listener can hear a realistic, 
stereophonic sound as if he were in the model space. 

[Second Embodiment] 

[0079] The sound environment simulator in accord- 
ance with the second embodiment is detailed in Fig. 13. 
The essential construction thereof is identical with the 
simulator of the first embodiment except that the sound 
environment analyzing unit 1 includes the echotime pat- 
tern computing unit 1 5 and impulse response computing 



unit 21; that the impulse response data for every location 
are computed with a computing server and stored into 
a memory in advance; and that the sound field repro- 
ducing unit 2 reads out the impulse response data and 

5 a head transfer function corresponding the location data 
from the memory to reproduce the stereophonic data 
from the dry source on real time. Otherwise, like com- 
ponents are labeled with like reference numerals with 
respect to the first embodiment, and the description of 

10 these component is not repeated. 

[Third Embodiment] 

[0080] The sound environment simulator in accord- 
?5 ance with the third embodiment is detailed in Fig. 14. 
The essential construction thereof is identical with the 
simulator of the first embodiment except that the sound 
environment analyzing unit 1 includes the head transfer 
storing unit 23 and a composite response data comput- 
20 ing unit 26. The composite response data computing 
unit 26 convolutes the head transfer function data and 
impulse response data to compute the composite re- 
sponse data in advance, and accordingly, the sound re- 
producing unit 24 convolutes the composite response 
25 data for each head direction into the dry source on real 
time with the DSPs. Otherwise, like components are la- 
beled with like reference numerals with respect to the 
first embodiment, and the description of these compo- 
nent is not repeated. 
30 [0081] With the simulator constructed as above, the 
number of the DSPs can be decreased as the number 
of incident directions increases. At the same time, less 
the data amount to be transferred to the DSPs becomes, 
faster the sound environment simulator can correspond 
35 to the dynamic changes of the listener. 

[0082] As has been explained, computing the time se- 
ries data of the arrival volume enables a precise predic- 
tion of the reverberant time. It also allows an accurate 
computation of the frequency characteristics of a rever- 
se berant sound. Thus, one or a few DPSs can also serve 
as a delaying circuit when the time length of the rever- 
berant sound and the frequency characteristics are an- 
alyzed and approximated. As a result, the realistic 
sound field can be reproduced by using less DSPs with- 
45 out impairing the listener's feeling of his presence in the 
model space. 

[0083] The reflected sounds are divided into the lower 
order reflections and higher order reflections and ana- 
lyzed by the respective computation methods to en- 

50 hance the work efficiency of the memory. However, the 
algorithm used for the higher order reflections can be 
used to all of the reflected sounds. 
[0084] In all the three embodiments, the sound envi- 
ronment simulator is designed to simulate the sound en- 

55 vironment which adequately corresponds to the chang- 
es occurring when the listener sits, stands up, walks 
around or moves his head. However, it may be designed 
to apply to a specific situation alone. 
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[0085] If the sound environment simulator is to be ap- 
plied to a case where the listener do not change his po- 
sition, the sound receiving point R can be fixed to a sin- 
gle point, which obviates the location data detecting unit 
31 and location data input unit 22, contributing to sim- s 
plifying the construction. 

[0086] If the sound environment simulator is to be ap- 
plied to a case when the listener moves his head in the 
same position, the sound receiving point R can be fixed 
to a single point as well, and the location data detecting io 
unit 31 is only to detect the directional change of the 
listener's head. Hence the amount of data used for com- 
putation decreases, which results in enhancing the 
overall operation speed. 

[0087] If the sound environment simulator is to apply *5 
to a case when the listener changes his position as he 
walks around, the height from the floor level to the lis- 
tener's ears have a constant value. Therefore, the three 
dimensional location data can be replaced with two di- 
mensional data. Hence, the overall operation speed can 20 
be enhanced for the same reason described in the pre- 
vious paragraph. 

[0088] With the sound environment simulator of the 
present invention, there may be a plurality of virtual 
sound sources, and even when they change their loca- 25 
tions, the realistic, stereophonic sound can be repro- 
duced as well; for it reproduces an adequate sound field 
by adding the impulse response data for each incident 
direction accordingly. 

[0089] The confined space was created as the virtual 30 
space in all the embodiments; however, it can be an 
open space if equations for the open space are used to 
compute the echotime patterns. Also, the present inven- 
tion enables the computation of the arrival volume of the 
sound of the higher order reflections as has been de- 3S 
scribed. Further, it can transduce the echotime pattern 
into the impulse response data in accordance with the 
listener's arbitrary position and head direction. Conse- 
quently, the present invention reproduces the realistic 
sound field by means of the existing VR equipment, 40 
which improves efficiency in designing for an audio and 
living room. 

[0090] Although the present invention has been fully 
described by way of example with reference to the ac- 
companying drawings, it is to be noted that various 45 
changes and modification will be apparent to those 
skilled in the art. Therefore, unless otherwise such 
changes and modifications depart from the scope of the 
present invention, they should be construed as being 
included therein. so 



Claims 

1. A sound environment simulator comprising: ss 

echotime pattern computing means (15) for 
computing an echotime pattern of a space; 



impulse response computing means (21) for 
computing an impulse response using said 
echotime pattern; 

reproduced sound generating means (25) for 
generating reproduced sound data using said 
impulse response; and 

output means (3) for outputting a sound by 
transducing said reproduced sound data into 
the sound; 

said simulator being characterised in that: 

said echotime pattern computing means (15) is 
adapted to compute said echotime pattern by 
a geometric method based on the theory of 
acoustics with input data as to conditions of 
said space; 

said impulse response computing means (21) 
is adapted to compute a response for each fre- 
quency band from said echotime pattern and to 
compute a composite response to generate 
said impulse response; and 
said reproduced sound generating means (25) 
is adapted to generate the reproduced sound 
data by convolution of said impulse response 
on non-reverberate sound source data. 

2. The simulator of claim 1 , wherein said simulator fur- 
ther comprises analysis data setting means (1 1 -1 4) 
for setting analysis data by receiving analysis con- 
dition data including a shape and physical data of 
walls of said space, dividing said walls into a set of 
certain sections, and setting a virtual sound source 
and a virtual sound receiving point at respective ar- 
bitrary positions within said space; and 

said echotime pattern means (15) is adapted 
to receive said setting analysis data and includes: 

time series data computing means (16) for 
computing, based on said setting analysis data, 
time series data related to arrival volume of a 
sound emanated from said virtual sound source 
to said virtual sound receiving point; and 
data transducing means (17) for transducing 
said time series data related to the arrival vol- 
ume into said echotime pattern by interpolating 
certain data. 

3. The simulator of claim 2 further comprising: 

dry source input means (24) for inputting said 
non-reverberate sound source data; and 
head transfer function storing means (23) for 
storing data related to a head transfer function; 
wherein said reproduced sound generating 
means (25) is adapted to generate the repro- 
duced sound data by convoluting said impulse 
response received from said impulse response 
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computing means (21) and said head transfer 
function received from said head transfer func- 
tion storing means (23) on said non-reverber- 
ate sound source data received from said dry 
source input means (24). 



outputting a sound by transducing said repro- 
duced sound data in the sound; 

said method being characterised by the step 

of: 



4. The simulator of claim 3, further comprising: 

location data detecting means (31) for detect- 
ing data related to at least one of a listener's 
head location and a listener's head direction; 
wherein said reproduced sound generating 
means (25) is adapted to extract said impulse 
response and said head transfer function cor- 
responding to said detected location data to 
convolute them on said non-reverberate sound 
source data. 

5. The simulator of claim 4, wherein said location data 
detecting means (31 ) is adapted to detect time var- 
ying data. 

6. The simulator of claim 4, wherein said output means 
(32) comprises a headphone. 

7. The simulator of claim 1 , further comprising: 



detecting data related to a listener's location; 
wherein the step of computing an echotime pat- 
tern computes said echotime pattern for a vir- 
10 tual sound source placed arbitrarily in the space 

by a geometric method based on the theory of 
acoustics; 

the step of computing an impulse response 
computes said impulse response by convolut- 

is ing said echotime pattern for each frequency 

band by a bandpass filter to compute a re- 
sponse for each frequency band, and by com- 
puting a composite response to compute said 
impulse response; and 

20 the step of generating reproduced sound data 

generates said reproduced sound data by re- 
ceiving a dry source of the sound source, and 
by convoluting said impulse response and a 
head transfer function corresponding to the de- 

25 tected listener's location on said dry source of 

the sound source. 



head transfer function storing means (23) for 
storing data related to a head transfer function; 
and 30 
composite response generating means (26) for 
computing a composite response for each di- 
rection of a listener's head by convolution of 
said impulse response and said head transfer 
function, 35 
wherein said reproduced sound generating 
means (25) is adapted to generate reproduced 
sound data by convolution of said composite re- 
sponse on the non-reverberate sound source 
data. 40 



10. The method of claim 9, wherein in the step of de- 
tecting data related to the listener's location, data 
related to dynamic changes of the listener's position 
and the direction of his head are detected as the 
location data. 

11. The method of claim 9, wherein in the step of de- 
tecting data related to the listener's location, data 
related to directional changes of the listener's head 
are detected as the location data. 

12. The method of claim 9, further comprising the step 
of: 



8. The simulator of claim 2, wherein said time series 
data computing means (16) is adapted to compute 
form factors to compute the arrival volume of the 
sound to the sound receiving point by using said 45 
form factors and a sound absorption coefficient of 
each section included in said physical data of the 
walls. 

9. A method of analysing a sound space, comprising so 
the steps of: 

computing an echotime pattern of a space; 
computing an impulse response using said 
echotime pattern; 55 
generating reproduced sound data related to a 
reproduced sound based on said impulse re- 
sponse; and 



dividing walls of the space into a set of sections; 
setting analysis conditions including the shape 
of the space, and physical data of the walls, the 
set of sections, a virtual sound source and a 
virtual sound receiving point; 
wherein the step of calculating said echotime 
pattern includes the steps of: 
computing time series data related to arrival 
volume of a sound emanated from said sound 
source to said sound receiving point for each 
incident direction and frequency band; and 
computing said echotime pattern by interpolat- 
ing certain data to said time series data. 

13. The method of claim 12, wherein the step of com- 
puting time series data includes the steps of: 
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computing data related to the arrival volume of 
the sound emanated directly from said sound 
source to said sound receiving point; 
computing data related to the arrival volume of 
the sound arriving at said sound receiving point 
by reflecting for up to a predetermined number 
of times by using a sound absorption coefficient 
and form factors of each section included in 
said physical data of the walls; 
computing data related to the arrival volume of 
the sound arriving at said sound receiving point 
by reflecting for more than the predetermined 
number of times by using the sound absorption 
coefficient and the form factors of each section 
included in said physical data of the walls; and 
computing time series data by dividing said ar- 
rival volume data into a set of groups by a cer- 
tain time interval. 

1 4. The method of claim 1 3, wherein in the step of com- 
puting data related to the arrival volume of the 
sound arriving at said sound receiving point by re- 
flecting for more than a predetermined number of 
times, a new virtual sound source is set by integrat- 
ing said arrival volume data of each group, and the 
arrival volume of sounds emanated from said new 
virtual sound source to any of the sections are re- 
peatedly computed. 

Patentanspruche 

1. Schallfeldumgebungssimulator, der aufweist: 

eine Echolaufzeit-Muster-Berechnungseinrich- 
tung (15) zum Berechnen eines Echolaufzeit- 
Musters eines Raums; 



zu berechnen; 

daG die Impulsansprechverhalten-Berech- 
nungseinrichtung (21) dazu angepaGt ist, ein 
s Ansprechen auf jedes Frequenzband auf dem 

Echolaufzeit-Muster zu berechnen und ein 
Komposit-Ansprechverhalten zu berechnen, 
urn das Impulsansprechverhalten zu erzeugen; 
und 

10 

wobei die Schallwiedergabe-Erzeugungsein- 
richtung (25) so angepaGt ist, urn die reprodu- 
zierten Schalldaten durch Konvolution des Im- 
pulsansprechverhaltens auf den nicht wider- 
75 hallenden Schallquellendaten zu erzeugen. 

2. Simulator nach Anspruch 1, wobei der Simulator 
weiterhin eine Daten-Analyseeinstelleinrichtung 
(11-14) zum Einstellen von Analysedaten durch 

20 Aufnehmen von Analysezustandsdaten, die Form- 
und physikalische Daten der Wande des Raums 
umfassen, Aufteilen der Wande in einen Satz von 
bestimmten Abschnitten und Einstellen einer virtu- 
ellen Schallquelle und eines virtuellen Schallauf- 

25 nahmepunkts an jeweiligen, wahlweisen Positionen 
innerhalb des Raums aufweist; und 
wobei die Echolaufzeit-Muster-Einrichtung (15) so 
angepaGt ist, um die Etnstellanalysedaten aufzu- 
nehmen, und umfaGt: 

30 

eine Zeitfolgedaten-Berechnungseinrichtung 
(16) zum Berechnen, basierend auf den Ein- 
stellanalysedaten, von Zeiffolgedaten, die sich 
auf ein Ankunftsvolumen von Schall beziehen, 
35 der von der virtuellen Schallquelle zu dem vir- 

tuellen Schallaufnahmepunkt abgestrahlt wird; 
und 



eine Impulsansprechverhalten-Berechnungs- 
einrichtung (21) zum Berechnen eines Impuls- 40 
ansprechverhaltens unter Verwendung des 
Echolaufzeit-Musters; 

eine Schallwiedergabe-Erzeugungseinrich- 
tung (25) zum Erzeugen reproduzierter Schall- 45 
daten unter Verwendung des Impulsanspre- 
chens; und 

eine Ausgabeeinrichtung (3) zum Ausgeben 
von Schall durch Umwandeln der reproduzier- 50 
ten Schalldaten in den Schall; 

wobei der Simulator dadurch gekennzeichnet 
ist, daG die Echolaufzeit-Muster-Berechnungs- 
einrichtung (1 5) so angepaGt ist, um das Echo- 55 
laufzeit-Muster durch ein geometrisches Ver- 
fahren basierend auf der Akkustik-Theorie mit 
Eingabedaten als Bedingungen fur den Raum 



eine Datenumwandlungseinrichtung (17) zum 
Umwandeln der Zeiffolgedaten, die sich auf 
das Ankunftsvolumen beziehen, in das Echo- 
laufzeit-Muster durch Interpolieren bestimmter 
Daten. 

I. Simulator nach Anspruch 2, der weiterhin aufweist: 

eine Trockenquelleneingabeeinrichtung (24) 
zum Eingeben der nicht widerhallenden Schall- 
quellendaten; und 

eine Kopfubertragungsfunktionsspeicherein- 
richtung (23) zum Speichern von Daten, die 
sich auf eine KopfObertragungsfunktion bezie- 
hen; 

wobei die Schallwiedergabe-Erzeugungsein- 
richtung (25) so angepaGt ist, um die reprodu- 
zierten Schalldaten durch Konvolution des Im- 
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pulsansprechens, das von der Impulsan- 
sprechverhalten-Berechnungseinrichtung (21 ) 
erhalten ist, undderKopfubertragungsfunktion, 
die von der Kopfubertragungsfunktions- 
speichereinrichtung (23) erhalten ist, auf den 5 
nicht widerhallenden Schalldaten, die von der 
Trockenquelleneingabeeinrichtung (24) aufge- 
nommen sind, zu erzeugen. 

4. Simulator nach Anspruch 3, der weiterhin aufweist: 10 

eine Lagedatenerfassungseinrichtung (31) 
zum Erfassen von Daten, die sich auf minde- 
stens eine Kopfstelle eines Horers und eine 
Kopfrichtung eines Horers beziehen; is 

wobei die Schallwiedergabe-Erzeugungsein- 
richtung (25) dazu angepaGt ist, das Impulsan- 
sprechen und die Kopfubertragungsfunktion 
entsprechend den erfaGten Lagedaten zu ex- 20 
trahieren, um sie auf die nicht widerhallenden 
Schalldaten zu konvolutieren. 

5. Simulator nach Anspruch 4, wobei die Lagedaten- 
erfassungseinrichtung (31 ) so angepaGt ist, um Da- 2s 
ten zu erfassen, die sich uber die Zeit variieren. 

6. Simulator nach Anspruch 4, wobei die Ausgabeein- 
richtung (32) einen Kopfhorer aufweist. 

30 

7. Simulator nach Anspruch 1 , der weiterhin aufweist: 

eine Kopfubertragungsfunktionsspeicherein- 
richtung (23) zum Speichern von Daten, die 
sich auf eine Kopfubertragungsfunktion bezie- 35 
hen; und 

eine Komposit-Ansprechverhalten-Erzeu- 
gungseinrichtung (26) zum Berechnen eines 
Komposit-Ansprechverhaltens fur jede Rich- 40 
tung des Kopfes eines Horers durch Konvolut- 
ion des Impulsansprechverhaltens und der 
' Kopfubertragungsfunktion, 

wobei die Schallwiedergabe-Erzeugungsein- *5 
richtung (25) dazu angepaGt ist, reproduzierte 
Schalldaten durch Konvolution des Komposit- 
Ansprechverhaltens auf den nicht widerhallen- 
den Schallquellendaten zu erzeugen. 

so 

8. Simulator nach Anspruch 2, wobei die Zeiffolgeda- 
ten-Berechnungseinrichtung (16) so angepaGt ist, 
um Formfaktoren zu berechnen, um das Ankunfts- 
volumen des Schalls an dem Schallaufnahmepunkt 
unter Verwendung der Formfaktoren und eines 55 
Schallabsorptionskoeffizienten jedes Abschnitts. 

der in den physikalischen Daten der Wande enthal- 
ten ist, zu berechnen. 



9, Verfahren zum Analysieren eines Schallraums, das 
die Schritte aufweist: 

Berechnen eines Echolaufzeit-Musters eines 
Raums; 

Berechnen eines Impulsansprechverhaltens 
unter Verwendung des Echolaufzeit-Musters; 

Erzeugen von Schallwiedergabedaten, die zu 
einem reproduzierten Schall basierend auf 
dem Impulsansprechverhalten in Bezug ge- 
setzt sind; und 

Ausgeben von Schall durch Wandeln der repro- 
duzierten Schalldaten in den Schall; 

wobei das Verfahren gekennzeichnet ist 
durch die Schritte: 

Erfassen von Daten, die sich auf die Stelle ei- 
nes Horers beziehen; 

wobei der Schritt eines Berechnens eines 
Echolaufzeit-Musters das Echolaufzeit-Muster 
aus einervirtuellen Schallquelle, die wahlweise 
in dem Raum plaziert ist, durch ein geometri- 
sches Verfahren, basierend auf der Akustik- 
Theorie, berechnet; 

wobei der Schritt eines Berechnens eines Im- 
pulsansprechens das Impulsansprechen durch 
Konvolutieren des Echolaufzeit-Musters fur je- 
des Frequenzband durch einen BandpaGfilter 
berechnet, um ein Ansprechverhalten fur jedes 
Frequenzband zu berechnen, und durch Be- 
rechnen eines Komposit-Ansprechverhaltens, 
um das Impulsansprechverhalten zu berech- 
nen; und wobei 

der Schritt der Erzeugung reproduzierter 
Schalldaten die reproduzierten Schalldaten 
durch Aufnehmen einer Trockenquelle der 
Schallquelle erzeugt, und durch Konvolutieren 
des Impulsansprechverhaltens und einer Kopf- 
ubertragungsfunktion entsprechend der 
erfaGten Stelfe des Horers auf der Trocken- 
quelle der Schallquelle erzeugt. 

10. Verfahren nach Anspruch 9, wobei in dem Schritt 
eines Erfassens von Daten, die sich auf die Stelle 
des Horers beziehen, Daten, die sich auf dynami- 
sche Anderungen in der Position des Horers und 
der Richtung seines Kopfs beziehen, als die Lage- 
daten erfaGt werden. 

11. Verfahren nach Anspruch 9, wobei in dem Schritt 
eines Erfassens von Daten, die sich auf die Stelle 
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des Horers beziehen, Daten, die sich auf richtungs- 
maGige Anderungen des Kopfs des Horers bezie- 
hen, als die Lagedaten erfaGt werden. 

Vertahren nach Anspruch 9, das weiterhin den 5 
Schritt aufweist: 

Untertetlen von Wanden des Raums in einen 
Satz von Abschnitten; 

10 

Einstellen von Analysebedingungen, die die 
Form des Raums, und physikalische Daten der 
Wande, den Satz von Abschnitten, eine vtrtuel- 
le Schallquelle und einen virtuellen Schallauf- 
nahmepunkt umfassen; is 



von Gruppen durch ein bestimmtes Zeitinter- 
vail. 

14. Vertahren nach Anspruch 13, wobei in den Schritt 
eines Berechnens von Daten, die sich auf das An- 
kunftsvolumen des Schalls beziehen, der an dem 
Schallaufnahmepunkt durch Reflektieren mehr als 
eine vorbestimmte Anzahl von Malen ankommt, ei- 
ne neue, virtuelle Schallquelle durch Integrieren der 
Ankunftsvolumendaten jeder Gruppe eingestellt 
wird, und wobei das Ankunftsvolumen des Schalls, 
der von der neuen, virtuellen Schallquelle zu ir- 
gendeinem der Abschnitte abgestrahlt wird, wieder- 
holt berechnet wird. 



wobei der Schritt eines Berechnens des 
Echolaufzeit-Musters die Schritte umfaf3t: 

Berechnen von Zeiffolgedaten, die sich auf ein M 
Ankunftsvolumen eines Schalls beziehen, der 
von der Schallquelle zu dem Schallaufnahme- 
punkt fur jede einfallende Richtung und jedes 
Frequenzband abgegeben ist; und 

25 

Berechnen des Echolaufzeit-Musters durch In- 
terpolieren bestimmter Daten zu den Zeiffolge- 
daten. 

Vertahren nach Anspruch 12, wobei der Schritt ei- 30 
nes Berechnens der Zeiffolgedaten die Schritte um- 
faGt: 

Berechnen von Daten, die sich auf das An- 
kunftsvolumen des Schalls beziehen, der direkt 35 
von der Schallquelle abgegeben ist, zu dem 
Schallaufnahmepunkt; 

Berechnen von Daten, die sich auf das An- 
kunftsvolumen des Schalls beziehen, der an 40 
dem Schallaufnahmepunkt ankommt, durch 
Reflektieren von bis zu einer vorbestimmten 
Anzahl von Malen unter Verwendung eines 
Schallabsorptionskoeffizienten und von Form- 
faktoren jedes Abschnitts, die in den physikali- *s 
schen Daten der Wande umfaGt sind; 

Berechnen von Daten, die sich auf das An- 
kunftsvolumen des Schalls beziehen, der an 
dem Schallaufnahmepunkt ankommt, durch so 
Reflektieren mehr als die vorbestimmte Anzahl 
von Malen unter Verwendung des Schallab- 
sorptionskoeffizienten und der Formfaktoren 
jedes Abschnitts, die in den physikalischen Da- 
ten der Wande umfaflt sind; und 55 

Berechnen von Zeiffolgedaten durch Untertei- 
len der Ankunftsvolumendaten in einem Satz 



Revendications 

1. Simulateur d'environnement sonore comprenant : 

des moyens de calcul de configuration de 
temps d'echo (15) pour calculer une configura- 
tion de temps d'echo d'un espace ; 
des moyens de calcul de reponse impulsionnel- 
le (21 ) pour calculer une reponse impulsionnel- 
le en utilisant ladite configuration de temps 
d'echo ; 

des moyens de generation de son reproduit 
(25) pour generer des donnees de son repro- 
duit en utilisant ladite reponse impulsionnelle ; 
et 

des moyens de sortie (3) pour sortir un son en 
convertissant lesdits donnees de son reproduit 
en son ; 

ledit simulateur etant caracterise en ce que : 

lesdits moyens de calcul de configuration de 
temps d'echo (15) sont adaptes pour calculer 
ladite configuration de temps d'echo par un pro- 
cede geometrique base sur la theorie de 
I'acoustique avec des donnees d'entree quant 
aux conditions dudit espace ; 
lesdits moyens de calcul de reponse impulsion- 
nelle (21) sont adaptes pour calculer une re- 
ponse pour chaque bande de frequences a par- 
tir de ladite configuration de temps d'echo et 
pour calculer une reponse composite pour ge- 
nerer ladite reponse impulsionnelle ; et 
lesdits moyens de generation de son reproduit 
(25) sont adaptes pour generer les donnees de 
son reproduit par la convolution de ladite repon- 
se impulsionnelle sur les donnees de source de 
son sans reverberation. 

2. Simulateur selon la revendication 1 , dans lequel le- 
dit simulateur comprend, de plus, des moyens 
d'etablissement de donnees d'analyse (11 a 14) 
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pour etablir des donnees d'analyse en recevant les 
donnees de condition d'analyse comprenant une 
forme et des donnees physiques des parois dudit 
espace, pour diviser lesdites parois en un ensemble 
de certaines sections, et pour etablir une source de 5 
son virtuelle et un point de reception du son virtuel 
aux positions arbitrages respectives dans ledit 
espace ; et lesdits moyens de configuration de 
temps d'echo (15) sont adaptes pour recevoir les- 
dites donnees d'analyse etablies et comprennent : 10 

des moyens de calcul de donnees serielles 
dans le temps (16) pour calculer, en se basant 
sur lesdites donnees d'analyse etablies, des 
donnees serielles dans le temps en rapport '5 
avec le volume a I'arrivee d'un son provenant 
de ladite source de son virtuelle audit point de 
reception du son virtuel ; et 
des moyens de conversion de donnees (17) 
pour convertir lesdites donnees serielles dans 20 
le temps en rapport avec le volume a I'arrivee 
en ladite configuration de temps d'echo en in- 
terpolant certaines donnees. 

Simulateur selon la revendication 2, comprenant, 25 
de plus : 

des moyens d'entree de source seche (24) 
pour entrer lesdites donnees de source de son 
sans reverberation ; et des moyens de memo- 30 
risation de fonction de transfer! de tete (23) 
pour memoriser des donnees en rapport avec 
une fonction de transfert de tete ; 
dans lequel lesdits moyens de generation de 
son reproduit (25) sont adaptes pour generer 35 
les donnees de son reproduit par la convolution 
de ladite reponse impulsionnelle recue en pro- 
venance desdits moyens de calcul de reponse 
impulsionnelle (21) et de ladite fonction de 
transfert de tete recue en provenance desdits 40 
moyens de memorisation de fonction de trans- 
fert de tete (23) sur lesdites donnees de source 
de son sans reverberation recues en provenan- 
ce desdits moyens d'entree de source seche 
(24). 45 

Simulateur selon la revendication 3, comprenant, 
de plus : 

des moyens de detection de donnees d'empla- 50 
cement (31 ) pour detecter des donnees en rap- 
port avec au moins I'un de I'emplacement de la 
tete d'un auditeur et de la direction de la tete 
d'un auditeur ; 

dans lequel lesdits moyens de generation de 55 
son reproduit (25) sont adaptes pour extraire 
ladite reponse impulsionnelle et ladite fonction 
de transfert de tete correspondant auxdites 



donnees d'emplacement detectees pour les 
convoluer sur lesdites donnees de source de 
son sans reverberation. 

5. Simulateur selon la revendication 4, dans lequel 
lesdits moyens de detection de donnees d'empla- 
cement (31) sont adaptes pour detecter des don- 
nees variant dans le temps. 

6. Simulateur selon la revendication 4, dans lequel 
lesdits moyens de sortie (32) comprennent un 
ecouteur. 

7. Simulateur selon la revendication 1 , comprenant, 
de plus : 

des moyens de memorisation de fonction de 
transfert de tete (23) pour memoriser des don- 
nees en rapport avec une fonction de transfert 
de tete ; et 

des moyens de generation de reponse compo- 
site (26) pour calculer une reponse composite 
pour chaque direction de la tete d'un auditeur 
par la convolution de ladite reponse impulsion- 
nelle et de ladite fonction de transfert de tete, 
dans lequel lesdits moyens de generation de 
son reproduit (25) sont adaptes pour generer 
des donnees de son reproduit par la convolu- 
tion de ladite reponse composite sur les don- 
nees de source de son sans reverberation. 

8. Simulateur selon la revendication 2, dans lequel 
lesdits moyens de calcul de donnees serielles dans 
le temps (16) sont adaptes pour calculer des fac- 
teurs de forme pour calculer le volume a I'arrivee 
du son au point de reception du son en utilisant les- 
dits facteurs de forme et un coefficient d'absorption 
du son de chaque section comprise dans lesdites 
donnees physiques des parois. 

9. Procede d'analyse d'un espace sonore, compre- 
nant les etapes consistant a : 

calculer une configuration de temps d'echo 
d'un espace ; 

calculer une reponse impulsionnelle en utilisant 
ladite configuration de temps d'echo ; 
generer des donnees de son reproduit en rap- 
port avec un son reproduit sur la base de ladite 
reponse impulsionnelle ; et 
sortir un son en convertissant lesdits donnees 
de son reproduit dans le son ; 

ledit procede etant caracterise par I'etape consis- 
tant a : 

detecter des donnees en rapport avec I'empla- 
cement d'un auditeur ; 
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dans lequel Pdtape consistant a calculer une 
configuration de temps d'echo calcule ladite 
configuration de temps d'echo pour une source 
de son virtuelle placee arbitrairement dans I'es- 
pace par un procede geometrique base sur la s 
theorie de i'acoustique ; 
I'etape consistant a calculer une reponse im- 
pulsionnelle calcule ladite reponse impulsion- 
nelle par la convolution de ladite configuration 
de temps d'echo pour chaque bande de fre- 10 
quences par un filtre passe-bande pour calcu- 
ler une reponse pour chaque bande de fre- 
quences, et par le calcul d'une reponse com- 
posite pour calculer ladite reponse 
impulsionnelle ; et 15 
Petepe consistant a g6nerer des donnees de 
son reproduit genere lesdits donnees de son 
reproduit par la reception d'une source seche 
de la source de son, et par la convolution de 
ladite reponse impulsionnelle et d'une fonction 20 
de transfert de tete correspondant a Pemplace- 
ment detecte de I'auditeur sur ladite source se- 
che de la source de son. 

1 0. Procede selon la revendication 9, dans lequel, dans 25 
I'etape consistant a detecter des donnees en rap- 
port avec I'emplacement de I'auditeur les donnees 

en rapport avec les changements dynamiques de 
la position de I'auditeur et de la direction de sa tete 
sont detectees comme les donnees d'emplace- 30 
ment. 

11. Procede selon la revendication 9, dans lequel dans 
I'etape consistant a detecter des donnees en rap- 
port avec I'emplacement de I'auditeur, les donnees 35 
en rapport avec les changements de direction de la 
tete de I'auditeur sont detectees comme les don- 
nees d'emplacement. 

12. Procede selon la revendication 9, comprenant, de <o 
plus, I'etape consistant a : 



en interpolant certaines donnees a partir des- 
dites donnees serielles dans le temps. 

13. Procede selon la revendication 12, dans lequel 
I'etape consistant a calculer des donnees serielles 
dans le temps comprend les etapes consistant a : 

calculer des donnees en rapport avec le volu- 
me a Parrivee du son provenant directement de 
ladite source de son audit point de reception du 
son ; 

calculer des donnees en rapport avec le volu- 
me a Parrivee du son arrivant audit point de re- 
ception du son apres un nombre predetermine 
de reflexions en utilisant un coefficient d'ab- 
sorption du son et les facteurs de forme de cha- 
que section comprise dans lesdites donnees 
physiques des parois ; 

calculer des donnees en rapport avec le volu- 
me a Parrivee du son arrivant audit point de re- 
ception du son apres un nombre de reflexions 
superieur au nombre predetermine en utilisant 
le coefficient d'absorption du son et les facteurs 
de forme de chaque section comprise dans les- 
dites donnees physiques des parois ; et 
calculer des donnees s6rielles dans le temps 
en divisant lesdites donnees de volume a Par- 
rivee en un ensemble de groupes par un certain 
intervalle de temps. 

14. Procede selon la revendication 13, dans lequel, 
dans I'etape consistant a calculer des donnees en 
rapport avec le volume a Parrivee du son arrivant 
audit point de reception du son apres un nombre de 
reflexions superieur au nombre predetermine, une 
nouvelle source de son virtuelle est etablie en inte- 
grant lesdites donnees de volume a Parrivee de cha- 
que groupe, et les volumes a Parrivee des sons pro- 
venant de ladite nouvelle source de son virtuelle a 
Pune quelconque des sections sont calcules a plu- 
sieurs reprises. 



diviser les parois de Pespace en un ensemble 
de sections ; 

etablir des conditions d'analyse comprenant la 4 $ 
forme de Pespace, et les donnees physiques 
des parois, Pensemblede sections, une source 
de son virtuelle et un point de reception du son 
virtue! ; 

dans lequel I'etepe consistant a calculer ladite so 
configuration de temps d'6cho comprend les 
etapes consistant a : 

calculer des donnees serielles dans le temps 
en rapport avec le volume a Parrivee d'un son 
provenant de ladite source de son audit point ss 
de reception du son pour chaque direction in- 
cidents et chaque bande de frequences ; et 
calculer ladite configuration de temps d'6cho 
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Fig. 11 
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